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Abstract 

We study eA scattering in a model where the total nuclear structure func- 
tion is a convolution of components Wn for an isolated nucleon and W^/a 
for a nucleus composed of point-nucleons. Wn is represented by its elastic 
and lowest-energy inelastic parts, while Wn/a i s computed from its asymp- 
totic limit, supplemented by final state interactions due to binary collisions 
between the knocked-out and core nucleons. (q/M)W^^A P er nucleon appears 
to be practically independent of mass number A and momentum transfer q. 
Consequently, predicted inclusive cross sections per nucleon for given kine- 
matic conditions hardly depend on the target and the same can be extracted 
from the data. Exceptions are low energy loss regions, where relatively small 
cross sections are shown to be sensitive to details of the single nucleon mo- 
mentum distribution and to FSI, involving multiple collisions. The agreement 
with the data is good. 
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Over the past few years a substantial body of data has been taken on cross sections for 
inclusive scattering of electrons from various nuclei. We concentrate on the SLAC-Virginia 
data for the following incident electron energies and scattering angles [1-3]: e= 2.0 GeV, 



9 = 15°, 20°; e= 3.6 GeV, 9 = 16°, 20°, 25°, 30°; e= 4.0 GeV, 9 = 30°, each for a range of 



energy losses v. Those sets span 3- and 4-momentum transfers q 2 < 4.8 GeV 2 and Q 2 < 3.35 
GeV 2 respectively. Targets were He,C, Al, Fe and Au, while quasi-data for nuclear matter 
(NM) have been generated by extrapolation of the above to A — > oo [4]. Below we shall give 
a description, addressing in particular the physical content of this large body of data. 

Until now virtually all approaches are based on the Plane Wave Impulse Approxima- 
tion (PWIA), where for high q the interaction of the core with the knocked-out proton is 
neglected. Central in the PWIA is the single-hole spectral function. Except for NM [5], 
accurate results exist only for light nuclei (cf. [6]) and for general A one has to resort to 
approximations. The latter range from selected, parametrized one-hole and two-hole, single 
particle contributions [7], to a local density approximation based on NM results [8]. Beyond 
the PWIA one eventually adds Final State Interactions (FSI) leading to a Distorted Wave 
Impulse Approximation (DWIA). That result, corrected for estimated colour transparency 
effects, has been compared with the above data [8-10]. 

In an alternative approach for NM we start from the large-g limit of the underlying 
structure function or response [11]. Contrary to the IA, there exists for non-relativistic 
(NR) dynamics, a systematic way to incorporate FSI contributions beyond the asymptotic 
limit [12]. Since the final expressions for responses contain only observables with no explicit 
reference to the underlying NR dynamics, those have been assumed to be valid in general. 

The above-mentioned approach for NM may be generalized for finite- A nuclei. Yet we 
prefer a somewhat different variant which we now expound. We start with the Rosenbluth 
cross section for inclusive scattering unpolarized electrons on a target of A composite nucle- 
ons 
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^) =4aV 2 cos 2 (tf/2)/Q 4 , (1) 
all J m 

where e', are the measured electron final energy and scattering angular volume. Wa are 
total nuclear structure functions, which have to be computed for a nucleus consisting of, in 
principle, composite nucleons. Those functions we represent below by a convolution of two 
parts: Wn/a for a nucleus made up of point-nucleons and Wn for a composite nucleon, each 
relating to different degrees of freedom. With £q the nucleon recoil energy we thus write 

a = 1,2) 



W itA (q, v) = J du'W N/A (q, + e q )W itN (q, v 1 ) (2) 

Expressions of the form (2) have been proven in special cases, for instance in the PWIA [13]. 
It is also a non-perturbative exact result for a NR quark-cluster model and the structure 
functions there refer to on-shell targets at rest. An adapted form has been conjectured to 
hold generally [14]. 

Eq. (2) permits the use of actual data for Wn, but for the kinematic conditions of the 
SLAC experiment it suffices to select those components where the nucleon remains in its 
ground state or is at most excited to a A. Then with i] = Q 2 /AM 2 and Ge,m, the standard 
elastic form factors one has 

W N = W^ + W$ el (3a) 

W$ = W^6(v - Q 2 /2M) (3b) 
GKQ 2 ) + vG 2 M (Q 2 ) 



1+7/ 



+ 2 v G 2 M (Q 2 )tg 2 (e/2) 



(3c) 



In contradistinction to Wn, I^n/a nas to be computed. We do so with emphasis on FSI 
due to collisions between the high-momentum knocked-out nucleon and low-momentum core 
nucleons. For NR dynamics Gersch et al. derived for the reduced response per nucleon [12] 



^,y) = (,/M)^ N/A = /^n(p)i?( ( /,-[,-^-^]) (4a) 

dp 

—^-n(p)R(q,y-p z ) (4b) 
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The above describes the reduced response <p(q } y) of a a nucleoli with momentum p which, 
after absorption of the 4-momentum q of the virtual photon, undergoes FSI contained in 
R. The expression is then folded into the distribution n(p) of the initial momenta of the 
nucleon. Eq. (4b) makes explicit the replacement of the energy transfer v by a (relativistic) 
scaling variable [15] 

y ~ {yw) ~ q r 2M + M )~ qV 2MJ> (5) 
(A) is some average separation energy, which is negligible except in the immediate neigh- 
bourhood of the Quasi-Elastic Peak (QEP) y ^ , v Q 2 /2M. 

FSI effects vanish in the asymptotic limit q — > oo and corresponds to R(q } y — p z ) — > 
6(y — p z ) [12] and therefore 

/•oo 

<p as ( y ) = 2tt / dppn(p) (6) 
■'\y\ 

Our major concern is the FSI factor R in (4) for finite q. For NR dynamics with interactions 
V between constituents there exists a systematic expansion in 1 /q of the dominant incoherent 
part of W N/A [12]. With v q = q/M 

t NR (q,y)= I £-e<y>i NR ( q ,s) (7a) 



L NR 



'<?, s) = dr t ri - sq) + — / dr 2 p 2 (r 1} r 2 ; r t - sq } r 

J l V q J 

j ds'V(ri — r 2 — s'q) — sV(ri — r 2 — sq)^j + .... 
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(7b) 

dr 1 p 1 (r 1 ,r 1 - sq)exp[tt(q, r-i, s)] (7c) 

The above expression contains density-matrices p n {fi^ f'\] f j) of any order n which are 
diagonal in all coordinates 2 < j < n, except in the coordinate 1 of the struck particle. For 
pi we use the Negele-Vautherin Ansatz [16]. 

Pl (r 1} r[) « p(SMs) « p( ri ) j j^- 3 e-V- s n(p) } (8) 

where S = (ri + r' 1 )/2, s = ri — r^; p(r) = p\(r } r) and n(p) is the single-particle momentum 
distribution. For n=2 



P2(ri,r 2 , r' l7 r 2 ) w p(r 2 ) Pl (r 1} r[)( 2 (r 1} r[; r 2 ), (9) 
with the 'off-diagonal' pair-distribution function ( 2 written as [12] 



(2(^1, r[; r 2 ) w ^g(r 1 , r 2 )g(r[,r 2 ) (10) 

The expression (10) modifies through ( 2 the Hartree approximation for p 2} but neglects 
the corresponding Fock part which has a longer range. 

Effects due to a strong repulsion in Vjvjv are as usual accounted for by a partial summation 
involving the following off-shell eikonal phase 



X(?,ri,s) = -if 



/ ds'V(r — s'q) — sV(r — sq) 
Jo 



:iv. 



For sufficiently high q, collisions of low multiplicity suffice and in the calculations to be 
reported, we limited ourselves to the simplest, binary collisions. In a first cumulant repre- 
sentation in Eq. (7c) reads [11] [r = r\ — r 2 )) 

Cl(q, n,s) = J drp{r 1 - r)( 2 (r, s)t(q, r^s) (12a) 
f = e «'(xi+X2) _ 1 ( 12 b) 

The phases Xj above correspond to the two terms in (11). xi there is an integral of V over 
a finite path segment and can be shown to be the off-shell analog of the standard on-shell 
phase in eikonal theory [17]. The same holds for the corresponding profile Ti = e tx — 1. For 
short range interactions it is approximately related to the on-shell analog (r = 6, z) 

f i(r, s) w 9(s - z)9(z)T 1 (b) (13a) 
ri(6) w — (1 -iT)A(b) (13b) 



Eq. (13b) is a standard parametrization of the profile in terms of <r tot , the total NN cross 
section, r, the ratio of real to imaginary part of the elastic NN amplitude and Aib) its range. 
Although rigorously Xi = ~ s ~§^Xi there appears to be no comparable simple expression for 
T 2 = e* X2 — 1. We approximate 



f = e «'(xi+X2) _ 1 ~ (1 + z ^ 2 ) 



5 e 



(14) 



and thus for , Eq. (12) 



tt(q } r 1} s) « -(a< ot (l - zr g )/2^} j J 2 6A g (6) 

/ ebp(&i -b } z t - z)( 2 (b, z } s) - 3/5(61 - 6, zi, s)( 2 (&, 5, s) 
.Jo 



(15) 



Here and in the following, brackets imply an average over p } n parameters (p, n densities are 
taken to be the same), thus 



(a to \l -ir)/2j = A' 1 



Zaf(\ - zr pp )/2 + Nv%(l - ir pn )/2 



(16) 



Eqs. (7c), (8) and (15) provide the elements for a calculation of the structure function Wn/a 
for a nucleus composed of point-nucleons 



M 

W N/A (q,u) = — <f>(q,y) 



<t>{<i,y) 



dp 

ds 



n(\p\)R[q,p z - y(q,v] 



R (9,V)= I 7^z lVS I dr 1 p(r 1 )exp[n(q,r 1 ,s)} 



(17a) 
(17b) 
(17c) 



Eqs. (17) are easily seen to correspond to a local density generalization of the corresponding 
FSI expression for NM used in [11], but with the parametrization for T 2 there replaced by 
the one implicit in (14). Combination of (1), (2) and (17) produces the total inclusive cross 
section (ICS) per nucleon 

/ f d P 



d 2 <r eA (e',0,Q 2 ) fda(e',9,Q 2 ] 



dttde' 



dv' 



dn J M J "' J (2tt) 

-2/ 



n(\p\)R[ q,p z - y(q,v') 



W 2 , N (q, u-u' + S q ) + 2tg 2 (9/2)W 1 , N (q, u-u' + S q 



:is) 



One notices that the energy loss in the above nucleon structure functions Wn is an inte- 
gration variable, whereas v = e — e' in the Mott cross section (Eq. (1)) is fixed. The two 
therefore do not generally combine into the inclusive cross section for a nucleon. This though 
is the case for the nucleon-elastic part (i.e. the component of Wn where the nucleon remains 
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unexcited) and which is our main concern. With barred quantities defined as in (2) one finds 
(cf. Eqs. (17) 

/,/2-r ,\N,el I / r ft~ \ N,el\ 

In the expression ultimately to be compared with experiment one adds to the elastic part 
(19) a nucleon inelastic contribution [9,18] 

dnde' ~ [dnde'J + [dnde'J ^ 

As remarked, the above is in essence a NR result with relativistic kinematics. However, in 
the total FSI phase 0, Eq. (15), there is no explicit reference to the underlying NR potential 
model. We then postulate its validity in the relativistic regime. 

We first address the iV-elastic part (19) in (20) and thus cf>A = (<?/M)Wn/a- The input 
for its calculation is relatively simple: 

i) Single-nucleon densities p(r), taken to be the same for p,n [19] 

ii) Single-nucleon momentum distributions n(p) 

e -(vlvi? + f -Wk) 2 
n 2 (p) = n 2 ,gr{p) + n 2 ,ex{p) (21) 

The distribution ni, with pi = ^J^pp] p 2 = V^Pi', fi = 0.03; pf=250 MeV, is one, indepen- 
dent of A [20]. For C and Fe we tested in addition a two-component form n 2 , derived from 
a parallel decomposition of the spectral function [7]. Finally for NM we exploited the result 
computed in [5]. 

iii) For the pair-distribution function in the Ansatz (10) for ( 2 we used g(r) for NM [21] 

iv) pp and pn scattering parameters as in Table I of Ref. [11] for cr tot ,r, Aib) in (13). 
With quite similar pp and pn scattering parameters and identical single nucleon densities, 

target specificity of Wn/a is mainly expected from single particle densities and momentum 
distributions in Eq. (17). 
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We have computed <p(q } y) from Eqs. (17) and (15) for C, Al, Fe and Au (He has been 
disregarded, because of the neglect of core recoil). Figs. la,b show <p(q } y) for Au and NM for 
some representative values q (in GeV)= 0.50, 2.2, oo. Reduced responses for other targets 
practically overlap with (f>Au- An insert displays for all targets <Pa(<1 = 2, 2, y) the immediate 
y < neighbourhood of the QEP. 

All reduced responses for relatively low q manifest modest FSI effects as asymmetries 
around y=0 in the large \y\ wings and a small shift in the position of the QEP towards 
negative y. Those effects vanish rather slowly for increasing q as is typical for, non-confining 
interactions with strong short-range repulsion [22]. 1 

To the extent that cf>A{q,y) is close to </>^ s (j/), the target non-specificity of the former is 
readily traced to the A-independent momentum distribution rii(p) } Eq. (21) in the expression 
(6) for (t>°£. In order to test sensitivity to n(p), we exploit the availability for C and Fe of 
the two distributions in (21). Fig. lc shows a < 10% effect in the QEP and additional 
differences in the wings. 

Finer details in <f> occur for large negative y } where the response in the binary collision 
approximation falls below a percent of the peak value at the QEP and occasionally even os- 
cillates around 0. Although (p exact > 0, separate components of a response or approximations 
need not be non-negative. The occasional oscillations in the binary collision approximation 
for <f> for large |j/|, are thus direct proof of the local insufficiency of that approximation. Also 
without that incontrovertible evidence, small FSI contributions due to binary and multiple 



1 Differences between responses for any A and NM are characteristically larger than the same for 
any two finite nuclei. It reminds one that distributions n(p) for heavy nuclei are not 'close' to the 
one for NM. We also note that in spite of different statistics and interactions between constituents, 
and hence entirely different single-constituent momentum distributions, responses of systems, as 
different as are finite nuclei and liquid He, have remarkably similar shapes. The same has been 
observed by G. West [23]. 
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collisions are expected to compete for sufficiently large \y\. 

Next we consider ICS ratios per nucleon under identical kinematic conditions. Using (19) 
and the near-universality of the reduced responses cf>A for nuclei composed of point-nucleons, 
one has 

f th ( a = d2(J eA 1 /A 1 _ 

(l VeA 2 / A 2 <pA 2 

The same ratio can also be extracted from data, and for Ai, A 2 > 12 one finds confirmed that, 
within a 10-15 % margin, Ca^ A 2 ~ 1- This is a remarkable and vital piece of information 
because, contrary to the theoretical expression (19), the left-hand side of (22) depends in 
principle on all kinematic variables e, 6, v, and on the mass numbers. Disregarding scatter in 
data where smooth behaviour is expected, systematic deviations occur only for near-elastic, 
small energy transfer v regions, or equivalently (for fixed q) regions with large negative y. 
To those we shall return shortly. 

A special case is A 2 — > (N). The appropriate ratio (22) then resembles a scaling function 
and its nucleon-elastic part (y < 0) is from (19) seen to be 

tn,el _ d 2 (J e A/A M 

£a,N = , 2 -N,el = (23) 
d<T e{N) $ 

Figs. 2a, b show for Fe separately, and for C, Al, Fe and Au together the y < parts of the 
ratio ^a,jv(<?,J/) extracted from the 4 data sets with e =3.6 GeV. Those are seen to display 
reasonable quality, universal scaling, i.e. approximate independence on A as well as on q. 
Fig. 2c displays £ for NM; its relevance will be discussed below. For NM, as well as for real 
targets, the quality of scaling in y (5) appears to be appreciably better than in a scaling 
variable, appropriate to the PWIA [2,25]. 1 

The close agreement between experimental and computed ICS ratios does not imply the 
same for the cross sections themselves. In the study of the latter one is nevertheless guided 



1 Another special case is ^a,d-, Eq. (22), as function of Q 2 and the Bjorken variable x = Q 2 /2Mv 
for x > 1 (y ;$ 0) between the QEP and the elastic threshold. Particular emphasis has been placed 
on the range 1.9 > x > 1.4, where ^a,d reaches a plateau value 5-6, far in excess of ^A lt A 2 ~ 1 
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by the result (22), which strongly limits the information present in the complete body of 

(Icitci. 

Fig. 3a thus presents data for Al, Au, e = 3.6 GeV, 6 = 25°. Drawn and dashed curves 
are predictions for Al and Au respectively. Fig. 3b gives the ICS per nucleon for Fe, e = 3.6 
GeV, 6 = 16°. The solid and dashed curves there are results, using the two distributions 
in (21). One notes discrepancies in the region of small v which, for approximately fixed q } 
corresponds to the large negative y wing of Wn/a (cf. Figs. 1). That region has been shown 
to be sensitive to details in n(p) and to higher order FSI. The former is apparent in Fig. la 
and it is interesting that for the kinematic conditions for Fe relevant to Fig. 3b, the data 
fall in between the two predictions. 

Fig. 3c gives another set of Fe data for e = 3.6 GeV, 6 = 25°. Along our predictions 
(drawn and dashed curves as in Fig. 3b) we inserted the DWIA prediction of Sick et al [8]: 
For small v those overestimate the data by an order of magnitude, in contradistinction to 
our slight overestimate. The authors above report colour transparency, suppressing purely 
hadronic FSI (hatched there). 

Finally we show in Fig. 3d for e=3.6 GeV, 6 = 30° extrapolated pseudo-data as well as 
predictions for symmetric NM and observe good agreement down to the smallest measured 
v. In view of the above one is tempted to link this with the single-nucleon momentum 
distribution which has been accurately calculated for NM [5,28]. Here too we show the 
DWIA result [9] and the claimed colour transparency reduction. 



for any two targets with A\,Ai > 4. First, on purely kinematic grounds £a,d grows towards the 
D elastic limit x=2. Far more important is the off-the-mean spatial extension of D. Indeed the 
height of the plateau of £a,d simply reflects the dissimilarity for D and A of the density p and pi 
in (19): it contains no other information. Similar, though less outspoken deviations of ICS ratios, 
occur for ^,He referring to the uncharacteristically compact 4 He: ratios are ~ 1.5 for medium x, 
approaching ~2 for x 2 [3]. 
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The generally satislactory agreement between data and our predictions is as good as 
may be expected lor any calculation ol hadron observables with natural imperfections in the 
underlying theory as well as in the quality of input. 

We have already mentioned the PWIA as the conventional starting point for the treat- 
ment of inclusive scattering and one is drawn to a comparison. There appear to be clear 
advantages in the approach starting from the asymptotic limit. From a technical point of 
view, the input is far simpler, as is for instance manifest in the required momentum distribu- 
tion as opposed to the two-variable spectral functions. Of greater importance are essential 
differences: 

1) Our treatment is free of off-shell problems inherent to perturbative PWIA. 

2) Starting from the DWIA [8,9] no simple argument explains the observed approximate 
universality of the data. 

3) The relative contribution of FSI to the total ICS is appreciably larger for the PWIA 
than for the asymptotic limit as starting point (cf. Fig. 3 in [9]). This is also reflected 
in the quality of scaling. In using a relativistic West scaling variable j/, one apparently 
includes part of the FSI beyond the PWIA approximation (see [30] for a recent example in 
the relativistic domain). 

In the end predictions for the \ow-v region based on hadron dynamics appear to differ 
considerably. The overestimate of the retained hadronic DWIA is reported to be off-set by 
colour transparency [9] (see below). 

We summarize our study of inclusive scattering of electrons from nuclei in the few-GeV 
region, where the basic total nuclear structure function has been represented as a folding 
of structure functions Wn of a nucleon at rest, and Wn/a f° r a nucleus composed of point- 
particles. 

We focussed on the nucleon-elastic region of Wn, for which the eA inclusive cross section 
factors in the elastic eN cross section and the response Wn/a(<?? u) f° r a nucleus composed of 
point-nucleons. The reduced response <f>A(q,y) = (q/M ) Wn/a(<?, v ) appears to be a nearly 
universal function, independent of A and of the momentum transfers q, leaving only a 
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modest role lor FSI over and above the asymptotic limit. An immediate consequence is the 
approximate target independence ol inclusive cross sections per nucleon, which is also borne 
out by experiment. A second consequence is the allocation of the strong variation in the 
eA cross section to the same in eN. The Ansatz (2) provides a simple explanation for the 
observed regularities. Calculations produce throughout satisfactory agreement with the ICS 
data. 

Regarding details we checked the influence of different single-nucleon momentum dis- 
tributions and found relatively large effects for low v ICS. We also emphasized the good 
agreement obtained for nuclear matter ICS down to the small-;/ region and mentioned the 
possible relation with a computed momentum distribution, available for that substance. 

A second agent which in principle contributes to ICS for low-;/ are FSI due to multiple 
collisions between the knocked-on and core nucleons. Those we emphasize, are part of the 
description using Glauber theory and a reliable calculation is definitely required. Rather 
than questioning the validity of Glauber theory for low-;/ ICS [8], we blame discrepancies 
there to imperfections in the application of a theory, formulated in terms of hadronic degrees 
of freedom. In particular we do not find evidence for colour transparency effects [8,9], nor 
are such effects expected for the relatively low Q 2 involved. 

The authors thank S.A. Gurvitz for a critical reading of the ms. 
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Figure captions 

Fig. la. The reduced structure function <f>A(q,y) = (q/M ) Wn/a(<?, ^), Eq. (17), for Au. 
The drawn, dashed and dotted curves correspond to, g(GeV)=0.5, 2.2, oo. The insert gives 
the small y < parts of </>^(2.2, y) for a number of targets. 

Fig. lb. Same as Fig. la for nuclear matter (NM). 

Fig. lc. Same as Fig. la for Fe, using the momentum distributions n\ (drawn line) and 
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n 2 (dashed line) in (22). 

Fig. 2a. The cross section ratio £a,n, Eq. (22), as lunction ol y lor Fe, e=3.6 GeV and 
4 scattering angles, showing empirical y-scaling. 

Fig. 2b. Same as Fig. 2a lor all targets, displaying approximate universal scaling. 

Fig. 2c. Same as Fig. 2a for symmetric NM. 

Fig. 3a. ICS, Eq. (20), for Al (squares), Au (crosses), e=3.6 GeV, 6 = 25°. Drawn and 
dashed curves correspond to Al and Au. 

Fig. 3b. Same as Fig. 3a for Fe, 6 = 16°. Drawn and dashed lines are predictions for 
the two momentum distributions rai, resp. n 2 in (22). 

Fig. 3c. Same as Fig. 3b for 6 = 25°. The hatched area gives the reported reduction of 
the dot-dashed DWIA results due to colour transparency [8]. 

Fig. 3d. Same as Fig. 3c for NM for e = 3.6 GeV, 6 = 30° [9]. 
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